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Abstract. We provide a general product formula for the solution of nonau- 
tonomous abstract delay equations. After having shown the convergence we 
obtain estimates on the order of convergence for differentiable history func- 
tions. Finally, the theoretical results are demonstrated on some typical nu- 
merical examples. 



1. Introduction 

Operator splitting is a widely used time discretization method for the numerical 
solution of complicated equations. The importance and main applications of these 
procedures is described, for example, in the monographs by Farago and Havasi [7], 
Holden et al. ^ and Lubich [TU| . 

The present paper investigates a special operator splitting for a class of nonau- 
tonomous delay differential equations. This method, which can be applied to equa- 
tions with distributed delays very effectively, was first investigated in Csomos and 
Nickel [5] and in Batkai, Csomos and Nickel [2] in the autonomous case. Recall 
from Bellen and Zennaro [4] that delay equations with distributed delay, especially 
those where the delay term is not separated from zero, are particularly difficult to 
solve numerically. Nevertheless, as we shall see, splitting methods work quiet well 
even in the latter case. 

To motivate this approach, let us consider the following equation. 



u{t) = b{t)u{t) + J n{t,a)u{t + a)da, t>s, 

-1 

u{s) = a; e R, 
[uis + a) = fia), C7e[-1,0], 

where b G C^(M), fi e L°°(R x [-1,0]), and t >^ fJ.{t,a) e Cj(R) for ah a G [-1,0]. 
In this case the delay operator ^(t) is defined by 



:= J ^iit,a)gia)da 
-1 

for ah g G Li([-1,0]). 

Choosing a time step h G (0, 1], first we start with xq :— x and /o := /. Then we 
set 



XI :=e'^''(^)(xo + /i*(s)/o) 
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(NCP) 



and 

fe(''+-W-)(a;o + /i$(s)/o), o e 

In the next step we repeat this procedure and replace xq with xi, /o with /i, and 
s with s + /i. Hence, we obtain an iteration process where in the /c*^ step we have 

(1) ] , . ^ ('e(''+'^)''(«+"(fc-i))(.x,._i+/i$(s + ;i(fc-l))/fc_i),ae [-/i,0], 

The aim of the present paper is to explain the convergence of this procedure also 
for more general equations. We do this by introducing an abstract setup allowing 
us a general convergence result of the procedure. Further, for differentiable initial 
function (i.e., classical solutions) we also obtain estimates on the order of conver- 
gence. In the following, we summarize some basic facts on product formulae for 
abstract evolution equations. Then in Section 2 we rewrite the nonautonomous 
delay equation as an abstract evolution equation and prove the convergence of a 
general product formula. Section 3 is devoted to the investigation of the order of 
convergence, and in Section 4 we present numerical examples demonstrating the 
power of this approach. 

First, let us recall some general facts about splitting of nonautonomous equa- 
tions. Consider an evolution equation of the form 

f^«(i) = (A(i)+BW)w(t), t>s6R, 
\u{s) = yC^X. 

We always suppose that this equation (jNCPp is well-posed, i.e., there is an evolution 
family W (also called semi-dynamical system) solving it. For well-posedness of 
nonautonomous evolution equations we refer to the surveys in Nagel and Nickel 
[TI]. Fazy Hg or Schnaubelt in [5]. 

Now, a splitting formula, as below, is especially useful, if we are able to solve 
effectively the autonomous Cauchy problems 

f^u{t) = A{r)u{t) 

^«(t) = B{r)v(t) 

with appropriate initial conditions for every fixed r. This is usually the case, if 
the operators A(r) and B(r) are partial differential operators with time dependent 
coefficients, or time dependent multiplication operators, and this is particulary so 
for delay equations considered in this paper. 

The following general convergence result can be proved, see Batkai et al. [TJ Theo- 
rem 4.2]. 

Theorem 1.1. Suppose the following: 

a) The nonautonomous Cauchy problem corresponding to the operators {A{-)+B{-)) 
is well-posed. We denote the evolution family solving (jNCP|) by W. 

b) The operators A(t) and B{t) are generators of Co-semigroups of type {AI,uj) 
(M >1 and u eR), and 

{lu, oo) C piA{t))r]p{B{t)) forallteM. 

and 



sup 

sGR 



ln< 



p—n 
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c) The maps 

R{\,A{t))y, t^R{X,B{t))y 
are continuous for all \ > uj and y G X . 
Then one has the convergence 

n-l 

(2) W{t,s)y^ lim FT (e'^^(^+''^e'^^(^+''^)y 

1=0 

for all y G X, locally uniformly in s,t with s <t. 

2. Splitting for the delay equation 
Consider the abstract delay equation in the following form: 

dt ' 



' ^w(t) = A{t)u{t) + <i>{t)ut, t > s, 



(3) 



u{s) ^ X e X, s eR, 

Us ^ f eL\[-l,0];X) 



on the Banach space X, where A{t) generates a strongly continuous contraction 
semigroup on X and : L^([— 1,0];X) — > X is a bounded and linear operator 
depending continuously on the parameter t G M. The history function ut is 
defined by Ut{a) := u{t + a) for cr e [—1,0]. Note that point delays are excluded 
from this context, but distributed delays, even those that live up to 0, are contained 
in this setting. 

In order to rewrite ([3]) as an abstract Cauchy problem, we take the product space 
£ :~ X X L^([— 1, 0]; X) equipped with 1-sum norm, and the new unknown function 
as 

t ^ U{t) :^ f E £. 



Ut 

Then ^ can be written as an abstract Cauchy problem on the space £ in the 
following way: 



(4) 



i-^u{t) - g{t)u{t), t > s. 



\ W(s) = (pe£, 

where the operator g(t) is given by the matrix 

(6) sm^^r-^w *1> 

\ " dtj 

on the domain 

DiGit)) :- {(p G DiA{t)) X W'-\[~1,0];X) : /(O) = x} . 

As in Batkai and Piazzera [21 Corollary 3.5, Proposition 3.9] one can show that the 
delay equation ([3]) and the abstract Cauchy problem (|4]) are equivalent, i.e., they 
have the same solutions. More precisely, the first coordinate of the solution of ^ 
always solves ([3]), i.e., 

U{t) = TTlUit), 

where tti is the projection to the first coordinate in £. Due to this equivalence, the 
delay equation is well-posed if and only if the operator g(t) generates an evolution 
family on the space £. 

Since the delay operators are bounded, the delay equation ^ is well-posed, 
which follows form a much more general well-posedness result by Hadd, Rhandi 
and Schnaubelt [5J Proposition 3.5]. That is, there is an evolution family W such 
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that for fixed s and t > s the function u'^^^t) = 7riW(t, s)(p is a solution of ([3]) 
for (p G D{Q{s)). In particular, for < s < t we have 



where fulfills w(")(t)=.T+ / A{r)u^''\r)dr + I <i>{r)ui''^dr 



and uf^(r) = + r) for r G [-1, 0], t + r > s. 

Furthermore, we have the next relation 

u'f'^ (r) = f{r + t — s) for t + r < s. 

Now we make the main assumptions implying the convergence of the splitting pro- 
cedure. In the autonomous case, i.e., when A{t) ~ A and — $, the following 
was investigated in the papers by Csomos and Nickel [5] and Batkai, Csomos, and 
Nickel H]. 

Assumption 2.1. a) The operators A{s) generate the strongly continuous con- 
traction semigroups {V'-^'' {t))t>o on X for all s £ R. 

b) D{A{s)) =: D for all s £ R and the function s i— ?■ R{l,A{s)) is continuous. 

c) The delay operators $(s) : L^([-l, 0]; X) X are bounded for all s e R. 

d) The function s H> <I>(s)/ is bounded and continuous for every / £ L^([— 1, 0]; X). 

Let us now describe in detail the approximation procedure we will apply. We split 
the operator in Q as 

git)^Ait) + Bit), 

where the sub-operators have the form 
(6) '^'^ 

Since A{r) is a generator and $(r) is bounded, the operators A{r) and B{r) generate 
the strongly continuous semigroups (T''"'(i))t>o and {S'^^'' {t))t>o, respectively. It 
is shown in Batkai and Piazzera [31 Theorem 3.25] that 7" is given by 



' [ v}^' m) ' 



where {T{t))t>o is the nilpotent left shift semigroup defined by 

10, if (T G [— t, OJ 

for all / e Li([-l,0];X), and v}''^ is 



l/W(t + cr)a;, if creH,0], 
0, if ae[-l,-t) 



for all a; G X. Since $(r) is a bounded operator, B{r) is also bounded on £. 
Therefore, the semigroup S{r) generated by B{r) takes the form 
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where /, /, and I denote the identity operators on X, L^{[~1,0];X), and £, re- 
spectively. We then have the following general convergence result explaining the 
convergence of the procedure described in ([T]) . 

Theorem 2.2. Under Assumvtion \2. 1\ the solution of the abstract delay equation 
^ is given by the formula 

Proof. The convergence is a direct consequence of Theorem 11.11 applied to the gen- 
erators A{r) and B{r). As mentioned above, by Hadd, Rhandi and Schnaubelt [5] 
the Cauchy problem associated to the operator A{r) +B{t) is well-posed. Applying 
Theorem 11.11 we only have to check the stability assumption. By using the same 
arguments as the ones appearing in the proof of Csomos and Nickel [Sj Theorem 
4.2], we get that 

lir« W(p \\s = II + \\v}'^x + r(t)/||i < ||.t|i + ||/||i + tM 

<(l + t)(||x|| + ||/||i) = (l + t)||(p||i, 

hence 

Iirww(;)ll£<i + i- 

One also obtains 

115^(011 <l + i|l'!>WII- 
From this the stability can be obtained: 

||n( nyjG " '/ '■') IN + ^) 

with w = 1 + sup,gjj ||(f>(s)||. □ 

3. Order of convergence 

We now investigate the order of convergence of the splitting method from the 
previous section. To this end, we have to make further assumptions on the operators 
involved, and of course, some regularity on the initial data need to be assumed, too. 

Assumption 3.1. 

a) The operator A{s) is bounded and generates the strongly continuous contraction 
semigroup {V''^^t))t>o on X. 

b) The delay operators $(s) : L^([-l, 0]; X) — ^ X arc bounded for all s G K. 

c) The function s i— >■ A(s)x is bounded and locally Lipschitz continuous, i.e., for 
all To > there is Ljo > such that 

\\A{s)x ^ A{t)x\\ < LTjx\\\t - s\ 

for all \t\, \s\ < To. 

d) The function s i— >■ $(s)/ is bounded and locally Lipschitz continuous, i.e., for 
all To > there is Ltq > such that 

ms)f~mf\\<LTo\\f\\i\t~s\ 

for ah \t\, \s\ < Tq. 

These assumptions enable us to show the first order of convergence for classical 
solutions. 
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Theorem 3.2 (Local error estimate). Let Tq > be fixed. Then there is a 
constant C > such that 

||r(^)(M5(^)(M(;) -wGs + /.,.)(;) II < ch^{\\x\\ + \\f\\, + \\r\\,) 

holds for alike [0,1], sG [-TcTq] anrf / G W^-^hl, 0]; X), a; = /(O). 
Proof. Recall from the previous section that for < s < s + /i we have 

i^'Hs + h)^ 



/u(''>(s + h)\ 



s+h 
s+h 



s+h 



where fulfills u^'^s + h) = x + / A{r)u^''\r)dr + / $(r)4'^dr 



and 



u|f^(r) = u^^'(/i + r) for/i + r>s. 



From this it follows that u^^^ : [s, s + 1] — X is Lipschitz continuous with constant 
^(ll^^ll + ||/||i) with L dependent only on \\A\\oo and ||$||oo- Let us calculate the 
product 



(7) 



f V(^Hh)x + hV^^\hMs)f \ 



and compare the first component here with u('^^(s + h). We can write 
V^'\h)x + hV^'\h)<i>{s)f - u^'\s + h) 

s-\-h s-\-h 

= V^'^h)x + hV'-''\h)<^>{s)f -X- [ A(r)J^'(r)dr - / <i>{r)ui''Ur, 



and by writing out the scries expansion of V^'^^ (h) wc obtain 
V^''\h)x + /iX^(")(/i)$(s)/ - m(^'(s + h) 

s-\-h s-\-h 

= x + hA{s)x + h<i>{s)f + 0{h^)-x- J A{r)u^''\r)dr - J $(r)4'^dr 

5 S 

(8) 

s-\-h s-\-h 

= hA{s)x~ J A{r)u'-''^{r)dr + 'l>{s)f - J ^{r)ui'Mr + 0{h^), 

S S 

where 0{h'^) denotes a term bounded in norm by C-/i^(||a:|| + ||/||i) with a constant 
C that depends only on the bounds of ||A||oo and ||$||oo- We now can write 

s-\-h s-\-h 

hA{s)x- [ A{r)u^''\r)dr < [ \\A{s)x - A{r)u^''\r)\\dr 



s-\-h s-^h 

< [ \\A{s)x ~ A{r)x\\dr + [ \\A{r)x ^ A{r)u^''^ {r)\\dr 



s+h 



< {L'\\x\\ + \\A\\^L) / {r - s)dr ^ Oih^) 



SPLITTING FOR DELAY EQUATION 



7 



where L' is the Lipschitz constant of A on [— TojTo + !]• A very similar reasoning 
works for the other two terms in ([5]): 

s+h S'\-h 

t^{s)f- j $(r)4'^ < j ||3'(s)/-$(r)4'^||dr 

S S 

s-\-h s-\-h 

< / ||$(s)/-$(r)/||dr+ / ||$(r)/-$(r)4^)||dr 



s+h 



s+h 



<L'||/||i / (r-.)dr+||$|U / \\f{a)-ui^Ha)\\dadr 



s-\-h s—r 



#!l/l|i/^' + ll*l|oo J J \\f{^)-f{<y 



r — s)\\d(7dr 



s+h 



+ ll*lloo / / \\f{a)^u^'\a + r)\\d<jdr 



s s—r 
s-\-h s—r 



+ 11*11. 



||/(a)-/(a + r-s)||dadr 



+ ||$||oomax{||/||o,,||7.(^)|[,,,+i]||o,} 
Since / G W^'^([— 1, 0]; X), we can continue the estimation: 



s+h 



s+h 



< 



L' I 



h' + moc. / ||/'||i(r-s)dr 



+ ||a>||oomax{||/|U,h(^^|[.,s+i]||oo}f 
<C(||x|| + ||/||i + ||/'||i);.^ 
By summing up we obtain the estimate: 

\\V^'Hh)x + hV^''){h)^s)f - u^'^h + s)\\< Ci\\x\\ + ll/lli + \\f\\i)h\ 
Hence the assertion for the first coordinate in ([7]) is proved. 

Let us now turn our attention to the second coordinate of ([7]). For r € [~1,0] we 
have the following: 

Eh+r>0 

(l^i^'x + hvl'^mf + T{h)f - u[%) (r) 

= V^''\h + r){x + h>^{s)f)+{)-u^''\s + h + r), 

and if ft, + r < 

{Vt'\ + hvt'^^{s)f + T{h)f - u'^:l^) (r) 

= - + /(ft + r - s) - /(ft + r - s) = 0. 

We estimate the L^-norm of 
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By using what is proved in the above for the first coordinate, the pointwise estimate 
for the integrand, we obtain that 

+ r)x 



< J \\V'-'\h + r)x + ih + r)V^''\h + r)<i>{s)f -u^'^s + h + r^dr 

-h 



+ / \\rV^''\h + r)<i>{s)f\\dr 





< C{\\x\\ + ll/lli + II/'Ili) J (/, + rfdr + I|<I>IU||/||iy 

= c{\\x\\ + \\fh + \\rh)h\ 

The proof is hence complete. □ 

We now can prove the first order convergence of the sequential splitting. Passing 
from local error estimates to convergence is done by the standard trick of telescopic 
summation. 

Theorem 3.3. For every Tq > there is constant C > such that for all f G 
W^''^{[~l,0];X) and x = /(O) the inequality 

holds for all s G [—To, Tq], < G [s, s + Tq] and for all n G N, where h ~ . 

Proof. Take tiq G N so large that Tq/jiq < 1 holds. Fix n > no, s as in the 
assertion and set h := {t — s)/n. By telescopic summation we obtain 

n-l 

(9) n '^^"^"''^ {h) - w{t, s) 

n—1 n—1 

^ Yi r(^+-'"") {h)s^^+^^^ (h) ~l[w{s + {j + i)h,s + jh) (p 

n—1 / n—1 
k=0 \ j = k+l 

X [r^'^+'^f'') {h)S^'+'''''> (h) - W{s + {k + l)h, s + kh)] X 

k-1 \ 

X (j[W{s + {j + l)h,s+jh)){^^)\ = 

j=0 ) 

n—1 / n—1 

= E ( n r^'+'''')(/i)5(^+^''')(/i))x 

X [r^<'+^^) {K) - W(.s + (fc + \)h, s + kh)] X 
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X W(s + fc/i,/i)(pj. 

+ 



we have ||a:fe||, ||/fe|| < (^11(^11. Therefore can we conclude from Theorem 13.21 that 

T^^+^''Hh)S^'+''''Hh)^Wis+ik + l)Ks + kh)Y2)\\<Ch'{M^ 
holds for all fc = 0, . . . , n — 1. From this and from ([9]) it follows 

n-l 

7-(.+.7.)(;^)5(.+,7.)(^)(p _ W(t,.)(p|| < nCh'{\\x,\\ + IIMIi + ll/^lli), 

hence the assertion. □ 

4. Numerical examples 

In this section we present our numerical examples obtained by the numerical 
code which applies the scheme described in Theorem 12.21 and in ([l} . The program 
code we apply is a modification of the code appearing in Csomos and Nickel [5]. 
In order to check the convergence of the numerical scheme in the nonautonomous 
case, and compare the solutions in the autonomous and nonautonomous cases, we 
will investigate the following examples. 

Example 4.1. Let X = R, B~beM. and consider 



^^u{t) = bu{t) + J ^i{t, a)u{t + a)da, t > 0, 

-1 

u{0) =x eR, 
uo = ./eLi([-l,0];R), 
for some /x G L°°(R x [—1,0]). In this case the delay operator is defined by 



u 



for all g £ L^([— 1, 0]; R). Let us choose the initial values as x = 1 and f{a) = 1 — a 
for a G [—1,0], and b = —1. As a particular example we choose the following 
functions fi: 

a) fj,{t, ct) = 1 in the autonomous case, and 

b) n(t,a) = 1 — sint in the nonautonomous case, 
for i > and a e [-1,0]. 

Example 4.2. Let us consider X^R, B^b^R and 

^M(t) = +/i(t)M(t - 1), t>0, 

u{0) = a: G M, 
uo = /gL1([-1,0];R). 
The delay operator in this case is 

Mi)g(-i) 

for all g G W^'^([— 1, 0]; M). Let choose the initial values again as a; = 1 and 
/(ct) = 1 — ct for (t G [—1,0), and b = —1. As above, we consider the functions fj, 
again: 
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a) /i(t,CT) = 1 in the autonomous case, and 

b) iJ,{t,a) = 1 — sint in the nonautonomous case, 
for i > and a e [-1,0]. 

Convergence of the numerical scheme. In order to examine the convergence of the 
numerical scheme described in Theorem l2.21 we plot the values of split solution m„ 
obtained by the formula 

/V^s+P^Hh) \ // h'f{s+ph)\ 

yis+ph) j.^^)) [q J ){f) 



p=0 



for different values of time step h. The results for the nonautonomous case are 
shown on Figure [TJ One can see that the split solutions converges to the exact 
solution if h decreases. The convergence of this numerical scheme in the autonomous 
case has been already investigated in Csomos, Nickel Section 5.3]. 




Figure 1. Results on the convergence of the numerical scheme for the 
nonautonomous delay equation with delay functions as in Example 14.11 (left 
panel) and Example 14. 21 (right panel). 



Long-time behaviour (difference between autonomous and nonautonomous cases). 
The long-time behaviour of split solutions m„ of the autonomous and nonautono- 
mous delay equations is shown on Figure [5] for the delay functions in Examples 14.11 
(left panel) and 14.21 (right panel). It can be clearly seen that in the case of the 
nonautonomous equation the difference in the delay functions does not play any 
qualitative role, because the effect of the function fi (i.e. the sin wave) suppresses 
it. 





Figure 2. Long-time behaviour of split solutions of the autonomous and 
nonautonomous delay equations with delay functions shown in Example 14.11 
(left panel) and Example 14.21 (right panel) . 
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Difference between delay functions in Examples \4-1\ and \4-2\ On Figure[3]the effect 
of the different delay functions are sliown in the autonomous and nonautonomous 
cases, respectively. As we have already seen, the effect of the delay function is 
suppressed by the sin wave of function fi in the nonautonomous case. In the au- 
tonomous case, however, the (structure of the) delay function plays an impor- 
tant role. 




Figure 3. Effect of the different delay functions in the autonomous (left 
panel) and nonautonomous (right panel) cases. 
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